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We perform advanced radiation leakage microscopy of routing dielectric-loaded plasmonic waveg- 
uiding structures. By direct plane imaging and momentum-space spectroscopy, we analyze the 
energy transfer between coupled waveguides as a function of gap distance and reveal the momentum 
distribution of curved geometries. Specifically, we observed a clear degeneracy lift of the effec- 
tive indices for strongly interacting waveguides in agreement with coupled-mode theory. We use 
momentum-space representations to discuss the effect of curvature on dielectric- loaded waveguides. 
The experimental images are successfully reproduced by a numerical and an analytical model of the 
mode propagating in a curved plasmonic waveguide. 



INTRODUCTION 

Confinement and propagation of surface plasmons in 
a metal circuitry have received considerable interest for 
their capability to transport data with a large band- 
width in compact structures and devices. Among the 
different geometries capable of routing the flow of sur- 
face plasmon, dielectric-loaded surface plasmon polari- 
ton waveguides (DLSPPWs) [lj— [sj have recently emerged 
as a potential plasmonic architecture that can be inte- 
grated seamlessly with current silicon-on-insulator (SOI) 
photonic circuits 0, [Hj and can sustain transfer of infor- 
mation at high date rates [6]. A DLSPPW is made of 
a rectangular dielectric rib deposited on a metal film or 
strip [7] . The surface plasmon is confined in the dielectric 
layer and typical cross-sections required for an optimum 
confinement of the mode compare well with state-of-the- 
art SOI waveguides operating in the telecom bands. De- 
spite dramatically higher losses, the advantage of such 
plasmonic platform is that the optical index of dielectric 
material used to confined the mode can be externally con- 
trolled to realize active DLSPPW-based devices [8l4l3j. 

Understanding device performance such as transmis- 
sion loss and surface plasmon mode profile greatly con- 
tributed to the development of DLSPPWs. Near-field 
optical microsco py and far-field leakage radiation mi- 
croscopy (LRM) [14[ are instrumental for visualizing the 
confinement and propagation details of surface plasmons 
supported by this geometry 0, [TH, [l6| . For thin metal 
films, LRM is an especially useful characterization tool. 
It provides a diffraction-limited snapshot of the mode 
profile in the structure, and, by conoscopic imaging en- 
ables to extract the wave-vector distribution. Effective 
indices of the different modes and interactions devel- 
oping in a g iven structure can thus be readily deter- 
mined 

In this work, we performed an in-depth analysis 



of leakage radiation images obtained for two different 
DLSPPW-based routing structures: linear couplers and 
bent waveguides. By simultaneously imaging the con- 
jugated aperture and field planes of the microscope we 
unambiguously quantify the degeneracy lift occurring for 
strongly interacting DLSPPWs and directly visualize the 
symmetry of the coupled modes. Furthermore, we mea- 
sured the wave-vector distribution of 90°-curved DLSP- 
PWs and show its evolution with bend radius. The exper- 
imental images are compared to numerical and analytical 
calculations. 



EXPERIMENTAL SECTION 
DLSPPW fabrication 

The waveguides considered in this work are de- 
picted schematically in Figs. [T] (a) to (c). A 
Poly(methylmethacrylate) (PMMA) ridge is defined by 
electron-beam lithography onto a 65 nm-thick gold film 
evaporated on a glass substrate. The thickness of all 
DLSPPWs fabricated in this study is fixed at £=560 nm 
and the width at w=600 nm. For such dimensions, the 
DLSPPW structures shown in Fig. [T] are single-mode at 
telecom wavelength {ijj. Scanning electron micrograph 
of the routing elements are shown in Figs. [I] (d) to (f). 



Characterization setup 

Optical characterization is performed by using a home- 
made leakage radiation microscope (LRM). The principle 
of this method is only recalled here, more detailed can be 
found in the review by Drezet et a/.[2l|. The schematic 
view of our experimental set-up is illustrated in Fig. El 
An incident tunable laser beam (fixed at A=1510 nm 
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FIG. 1. (a) Schematic view of a basic DLSPPW with a width 
w of 600 nm and thickness £=560 nm. (b) Schematic view of 
two coupled waveguides separated by a gap d. (c) Schematic 
view of a curved waveguide with a bend radius R c . (d), (e) 
and (f) are scanning electron micrographs of typical devices 
corresponding to the configurations (a), (b), and (c), respec- 
tively. 
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FIG. 2. Schematic of the leakage radiation microscope used in 
this study. Launching of the plasmon mode in the waveguides 
is obtained by focusing a laser operating at A=1510 nm with 
a long working distance objective. Surface plasmon leakages 
are collected by an oil immersion objective and are imaged on 
two cameras placed at the conjugated image (field) plane and 
Fourier (aperture) plane, respectively. 



hereafter) is focused by a 100 x microscope objective 
on the extremity of a DLSPP waveguide. The sharp dis- 
continuity defined by the polymer structure acts a local 
scatterer and produces a spread of wave- vectors, some 
of which resonant with the surface plasmon modes sup- 
ported by the geometry. By controlling the incident po- 
larization parallel to the longitudinal axis of the waveg- 



uide [22j, the DLSPPW mode can thus be readily ex- 
cited. Leakage radiation microscopy (LRM) provides 
a far-field imaging technique to directly visualize sur- 
face plasmon prop agation and investigate its fundamen- 
tal properties [14j, [23j . This method relies on the collec- 
tion of radiation losses occurring in the waveguide during 
propagation. These losses are emitted in the substrate 
at an angle phase-matched with the in-plane wave- vector 
of the SPP mode. In our LRM, the plasmon radiation 
losses are collected by an oil-immersion objective with a 
numerical aperture (N.A.) of 1.49. A tube lens focuses 
the leakages in an image plane (IP) conjugated with the 
object plane where an InGaAs infrared camera is record- 
ing the two-dimensional intensity distribution. Images 
recorded at this plane provide direct information about 
the propagation of the surface plasmon mode developing 
in the DLSPPW. To complete the analysis, we have also 
used the Fourier transforming property of the objective 
lens to access the wave- vector distribution of the emitted 
light. The angular distribution of the rays radiated in 
the substrate and collected by the objective lens is trans- 
formed to a lateral distribution at the objective back focal 
plane. Quantitative measurement of the complex surface 
plasmon wave- vector, or equivalent ly its complex effec- 
tive index, consists at measuring the radial distance of 
the rays with respect to the optical axis. Access to mo- 
mentum space was done by inserting a beam splitter in 
the optical path and a set of Fourier transforming relay 
lens. The lenses are forming an image conjugated with 
back focal plane of the objective. This Fourier plane (FP) 
contains the two-dimensional wave- vector distribution of 
the leakage radiation and momentum-space spectroscopy 
can thus be readily performed. 



Calibrating leakage radiation images 

Figures [31(a) and (b) show images of the leakage radia- 
tion intensity of a single-mode DLSPPW recorded in the 
conjugated image and Fourier planes, respectively. The 
plasmon mode is launched at the bottom of the structure 
and propagates up the waveguide with an exponentially 
decaying intensity. The image of Fig. [31(b) is a direct mea- 
surement of the wave- vector content of the intensity dis- 
tribution shown in Fig. [31(a). For both FP and IP images, 
a calibration was performed prior to any data extraction. 
For calibrating IP images, we used the known waveguide 
length as a standard. With the system magnification 
one pixel represents ~0.6 jam. For FP images, the N.A. 
of the objective was used to calibrate k x and k y axes. 
The largest ring in Fig. [31(b) represents the N.A.=1A9 
specified by the manufacturer. The central disk is the 
numerical aperture of the xlOO excitation objective at 
N.A .=0.52. One pixel provides a AN. A. ~ 0.012. 

The DLSPPW mode propagating along this straight 
waveguide presented in Fig. [3] can be defined by two pa- 
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FIG. 3. (a) Intensity distribution recorded at the image plane 
of a surface plasmon mode propagating along a straight DL- 
SPPW. The excitation spot is readily visible at the lower por- 
tion of the waveguide, (b) Corresponding wave- vector distri- 
bution. The central disk represents to the numerical aperture 
of the illumination objective (0.52). The DLSPP mode is 
recognized as a bright line at constant k y /k x (arrow), (c) 
Exponential fit of the experimental decaying plasmon inten- 
sity along the waveguide leading to a L spp —42.2 /im. (d) 
Lorentzian fit of the FP plasmon mode signature centered at 
f3'/k = 1.262. The full with at half maximum An e // is in- 
versely proportional to L spp . 
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FIG. 4. (a) Schematic of a linear DLSPPW coupler, (b) to 
(g) are direct-space images of the intensity distribution in lin- 
ear DLSPPW couplers with edge-to-edge distance d equals to 
670 nm, 600 nm, 440 nm, 350 nm, 240 nm and 90 nm respec- 
tively, (h) to (m) are the corresponding wave- vector distribu- 
tions unambiguously demonstrating the degeneracy lift of the 
coupled waveguides for small values of d. The images were 
cropped to display only the lower part of the Fourier plane. 



from direct-space analysis. 



MOMENTUM-SPACE SPECTROSCOPY OF 
LINEAR DLSPPW COUPLERS 



rameters: its effective index n e ff and propagation length 
L spp . n e ff is expressed by the phase constant ft and 
reads n e // = ft /k Q where k Q = 2tt/X. The propaga- 
tion length L spp = (2/3") _1 where ft' is the attenua- 
tion constant of the plasmon mode. A complex propaga- 
tion constant is then evaluated from these two constants: 
(3 = ft + ift f . 

L spp can be directly extracted from direct-space image 
by fitting an exponential decay of the intensity along the 
waveguide / = Ioex.p(—y/L sp p). Here a L spp =42.2 /xm 
is determined from the fit to the experimental data il- 
lustrated in Fig. EJc). Experimentally, the FP image 
displayed in Fig. [3^b) contains more information than 
a direct space image since the real part and imaginary 
part of the effective index can be directly measured. The 
signature of the mode is represented by a single line 
at a constant n e //=fc 2/ /fc o =1.262±0.006. The intensity 
measured along k y /k Q at k x =0 is related to the surface 
plasmon through the following formula [8]: I(k x ,k y ) oc 
\Ho(k x )\ 2 /[(k y - ft) 2 + (l/2L spp ) 2 }. H (k x ) is the k x - 
Fourier transform of the guided magnetic field at the ob- 
jective focal point. The imaginary part of the effective 
index is also estimated precisely through a Lorentzian 
fit. The width An e // of the ft jk Q line is a measure 
of the losses experienced by the plasmon mode and 
is thus inversely proportional to its propagation length 
L spp [24]. We obtain from the FP image a L spp =43.3 fim 
in fairly good agreement with L spp =42.2 /im measured 



We now demonstrate the added-value of performing 
momentum-space leakage radiation spectroscopy of a 
classical integrated plasmonic routing device: a linear 
DLSPPW coupler (Tfl, [H-[27| . The linear coupler geom- 
etry illustrated in Fig. H{a) consists of two parallel and 
identical DLSPP waveguides separated by an edge-to- 
edge distance d varying from 670 nm to 90 nm. This ele- 
mentar y co nfiguration is well-known from coupled-mode 
theory [28[. When the gap distance d is reduced the de- 
generate modes propagating in the uncoupled waveguides 
are split into symmetric and antisymmetric modes. New 
propagation constants ft s and ft as , respectively are thus 
characterizing the symmetric mode and the antisymmet- 
ric mode, respectively, and they critically depend on d. 
A beating of these two modes can be observed in leakage 
radiation microscopy 0,113 where a mode propagating 
in one waveguide can be totally transferred to the second 
after a coupling distance L c [29[. Figures H(b)-(g) qual- 
itatively show the evolution of the beating pattern and 
the coupling distance L c for decreasing separation dis- 
tances d. More interesting are the corresponding wave- 
vector distributions depicted in the series of FP images 
in Figs. HJh) to (m). When the DLSPPWs are weakly 
coupled (<i=670 nm), the Fourier content of Fig. H(h) 
strongly resembles that of a single DLSPP mode already 
shown in Fig. [3fb). When d is reduced, a clear splitting 
of the modes is observed indicative of a strong interaction 
between the waveguides. The parity of the modes can be 
readily determined from e.g. Fig.HJm). The asymmetric 



mode has an odd parity with two maxima centered on 
each side of k y /k Q axis. 

Figure [5] illustrates the benefit of performing 
momentum-space spectroscopy described here. Fig- 
ures 03(a) and (b) are leakage radiation images recorded 
at the conjugated Fourier plane and image plane for a lin- 
ear coupler with d=440 nm, respectively. The effective 
indices of the symmetric n s e ^ and antisymmetric n^ff 
modes are evaluated by Lorentzian fits of crosscuts of the 
momentum distribution along the k y /k axis marked by 
the circles in Fig. 03(a). n^ff was evaluated at two differ- 
ent wave- vector positions with respect to the k y /k Q axis, 

labeled as ft-"// an( ^ n eff- The ^ s ^° ^ ne data are repre- 
sented in Figs. [5](c) and (d) for the antisymmetric mode 
leading to n a e s ff ~ 7^=1.260^0.006. The effective index 
of the symmetric mode is measured at n s e ^=\. 310±0. 006 
(Fig. He)). 

The coupling length L c can now be evaluated using the 
following relation [29| 
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This procedure was repeated for the different DLSPPW 
separation distances d. The extracted values of the split 
modes and coupling length L c are reported in Table HI 
L c was estimated taking into account an average value 
of and n^ff. For comparison purposes, L c was also 
determined by analyzing the beating pattern recorded 
in the image plane (Fig. 03(b)). A longitudinal cross- 
section of the leakage intensity taken along each DLSP- 
PWs is shown in Fig. 03(f) f° r <i=440 nm. The transfer of 
the mode between the two waveguides is obtained after 
L c ~16 /im. Within our pixel resolution, this value com- 
pares well with the L c determined by the analysis of the 
degeneracy lift of the effective indices by Eq. [2] The ad- 
vantage of a momentum-space spectroscopy is that, un- 



like near-field measurement 



the effective indices 



can be directly measured and mode symmetry visual- 
ized [H, [H| • L c inferred from direct plane analysis only 
provides the difference between the two propagation con- 
stants \j3' s — j3' as \ without discriminating the symmetric 
mode from the antisymmetric one. 

To confirm these experimental results, we numerically 
simulated linear DLSPPW couplers with the commer- 
cial Finite-Element mode solver COMSOL. The opti- 
cal index of the PMMA waveguide is Tirp rfiraa = 1.535 
and that of the gold layer is n go id = 0.536 + i9.5681 at 
A=1510 nm [30[. The transversal electric field distribu- 
tion of the asymmetric and symmetric mode are shown 
Figs. [6] (a) and (b) for <i=400 nm. Figure [61(c) provides 
a comparison of the evolution of the degeneracy lift with 



FIG. 5. (a) and (b) are the respective leakage radiation 
Fourier and image planes obtained from a linear DLSPPW 
coupler with <i=440 nm. (c) and (d) are Lorentzian fits of 
the asymmetric mode at the location marked by the circle 
in (a), (e) is a Lorentzian fit of the symmetric mode, (f) 
Longitudinal intensity cross sections taken along the two cou- 
pled waveguides in (b) showing the energy transfer from one 
DLSPPW to the other defining the coupling length L c . 
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TABLE I. Measured values of the effectives indices of the 
asymmetric and symmetric modes and estimated coupling 
length L c for linear DLSPPW couplers with different edge- 
to-edge separation d. L c values were estimated independently 
by momentum-space spectroscopy and direct-plane analysis. 



separation distance d. The measured data quantitatively 
reproduce the response of the simulated device. 



MOMENTUM-SPACE SPECTROSCOPY OF 
CURVED DLSPPWS 

Experimental images 

In this section, we investigate the momentum dis- 
tribution of another well-known basic routing element: 
90° curved waveguides. This DLSPPW structure has 
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FIG. 6. (a) and (b) are the electric field profiles of the asym- 
metric and symmetric modes propagating in two coupled DL- 
SPPWs separated by <i=400 nm. (c) Comparison between 
the experimental splitting and the calculated values for vari- 
ous coupling distances d. 




FIG. 7. (a) and (b) Leakage radiation image of a bent waveg- 
uide with R c = 19 /im and its corresponding wave- vector dis- 
tribution, respectively. The Fourier content of the curved 
section is appearing as an arc linking the two k x and k y lines, 
(c) and (d) Computed intensity of the mode propagating in 
the routing element and its corresponding wave- vector distri- 
bution, respectively. 

ing along the y and x-oriented waveguides, respectively. 
The line at /c cc =constant is the input waveguide. This 
Fourier plane shows additional signatures such as the il- 
lumination wave-vector span (central disk) and planar 
plasmon modes supported by the Au/air or Au/PMMA 
taper interfaces. Of particular interest here is the Fourier 
signature of the bent section of the waveguide visible as 
an arc of circle linking the two k x and k y lines. 



Fourier plane model 



been extensively studied by various groups and the ef- 
fect of bend radius on the overall losses is well under- 
stood [IBI, l3ll - [34j . By performing momentum-space spec- 
troscopy of the supported mode, we visualize and an- 
alyze the wave-vector content of the bend section for 
this routing device. We demonstrate the limitations of 
momentum-space spectroscopy to extract modal prop- 
erties of this elementary building block. The curved 
DLSPPW is composed of two L=30 jiva long straight 
waveguides with iu=600 nm. These two input and out- 
put waveguides are linked by a circular 90° bend sec- 
tion of radius R c . Bent DLSPPWs with R c ranging from 
5 fim to 19 fim were fabricated. The leakage intensity 
distribution image of a R c =19 /im curved waveguide is 
shown Fig. Ufa). The corresponding wave- vector distri- 
bution is depicted in Fig[7Jb). The two straight lines at 
k x and k y constant are related to the mode propagat- 



Instead of full numerical simulations, we propose in 
the following a simple analytical model that explains 
the main features of the measured Fourier plane im- 
ages. To this aim, we approximated the mode that 
propagates in the bend structure by a gaussian profile 
with a finite propagation length. The characteristics 
of the gaussian profile was defined by the experimen- 
tal data obtained from a straight waveguide; namely 
n e ff= 1.262 and L spp =42.2 /im. The mode width was 
fixed at il>o=500 nm according to the analysis of Holm- 
gaard and Bozhevolnyi [2o| . 

The magnetic field is written as follow: 

• input straight guide (— L<x<0, y > 0) 



H y (x,y) = H Q ex.p 



-{y-Rcf 



exp[z/3x] 



(3) 
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FIG. 8. (a) Experimental momentum-space image of a R c = 19 /im bend waveguide. The dashed lines indicate the position 
of the momentum profile in the following graphs, (b) (c) and (d) show the comparison of the calculated and experimental 
wave- vector distributions along the k x , k y and -7r/4 axis, respectively. 



H Q exp 



exp[in e / fk Q x] ex-p[-x/2L spp ] , 



circular portion (x > 0, y > 0) 

H r (r,9) = H exp[-L/2L spp ] 
-(r - Re) 2 



(4) 



exp 



ex.p[ik b R c 9] , 



Here ws corresponds to the length over which the field 
outside of the waveguide decays by 1/e. With W3 c± wo, 
i?z ~ 13 /im. When the radius is below Ri, bending losses 
are induced by a radial displacement of the mode profile 
with respect to the waveguide axis [33|. This displace- 
ment pushes the field outside the waveguide leading to a 
lower phase velocity and a modification of the effective 
index 13511 . 



where r = y ' x 1 + y 2 , 6 = arctan and k b is 

the complex wave- vector of the mode in the curved 
section. 

• output straight guide (x > 0, — L < y < 0) 

-(x-R c ) 2 ^ 



H x (x,y) = H exp 



ex.p[-L/2L spp ] 



(5) 



ex.p[-R c 7r/4L spp ] ex.p[in eff k y} ex.p[-y/2L spp ] 



Considering curvature loss as an additional exponential 
decay is a good approximation to quantify the total loss 
induced by the bend and estimate the transmission of 
the 90° waveguide. However, this approximation is not 
representative of the real shape of the field along the 
bend and consequently cannot be used to model Fourier 
images of the kind displayed in Fig. Efb). 



The momentum representation H(k x , k y ) was obtained 
by a Fourier transform 

/fi c +3xw rR c +3xwo 
/ H(x,y) (6) 

-L J-L 

exp [— i(k x x + k y y)] dxdy. 

where the integration window is truncated to limit the 
calculation on the mode extension. Outside this area, 
the mode profile vanishes so that its contribution to the 
Fourier transform is negligibly small. 

To simulate the curved DLSPPW, we maintained the 
complex propagation constant in the curved section equal 
to that of the straight waveguides (fc& = f3). This as- 
sumption remains valid for R c > Ri where Ri is the lim- 
iting radius where bend losses can be neglected. Ri have 
been recently numerically estimated for DLSPPW (33| 
confirming an analytical expression historically used for 
standard optical waveguides [HI 

* > (7) 



For R c > Ri ~13 /im, we show that the propa- 
gation and momentum-space representations simulated 
with these basic assumptions (Figs. [7fc) and (d)) repro- 
duces the experimental observations of Figs. [T^a) and 
(b). We compare the profile of the calculated Fourier 
transform (Fig. Efd)) to the experimental wave- vector 
distribution (Fig.[7fb)) by extracting momentum profiles 
at different radial coordinates indicated in Fig. [8fa) for 
R c =19 /im. Figures Hfb), (c) and (d) show the calcu- 
lated and experimental wave-vector distributions taken 
along the k x and k y axis and at -7r/4, respectively. The 
position and width of the calculated wave- vector profiles 
are in very good agreement with the measured data val- 
idating thus our preliminary assumptions. Some extra 
signatures are visible on the experimental cuts: a peak 
at 1.49 in Fig. Efb) corresponds to the numerical aper- 
ture of the collection objective, and the contribution of 
the surface plasmon excited in the nearby gold film is vis- 
ible at 1.08 in Fig. Efd). Noticeable also is the presence 
of Gibbs oscillations revealed by this momentum- space 
spectroscopy (arrows Figs. [8] (b) and (c)) . The origin of 
these oscillations in reciprocal space is discussed below. 
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Analytical development 



In order to provide analytical expressions and propose 
a simple physical understanding of the measured Fourier 
images, we simplified further the model. We consider 
each part of the waveguide (straight and bend) indepen- 
dently from each others. The straight parts are approx- 
imatively 30 /im long, a length smaller than the longi- 
tudinal decay of the supported mode (L spp =42 /im). If 
we consider only the field in the input straight waveg- 
uide along the x coordinate, the Fourier transform of the 
magnetic field written in Eq. [3] reads 



H y (k x ) = H J exp[z(/3 - k x )x] dx , (8) 

where we neglected the ^/-gaussian profile for the sake 
of clarity but this could be easily considered. Then the 
intensity in the Fourier plane along k x direction writes 



\Hy(k x )\ 2 = \H \ 



2e" /? " L cos((/?' - k x )L) + , 



-20"L 



(9) 

If L -> oo then \H y {k x )\ 2 - k x ) 2 + p" 2 ) and 

follows a Lorentzian profile, as expected. If the straight 
part length L is smaller than L sppi the resulting Fourier 
transform is not simply defined by a Lorentzian function. 
Because of the finite integration limit, extra oscillations 
are becoming visible (Gibbs oscillations) as illustrated in 
Figs. |Sfb) and (c). For long waveguide (as in Fig. [3ja) 
where L > L spp ), no oscillation of the wave- vector dis- 
tribution is observed. Figure [9] shows the impact of the 
integration boundaries L on the momentum space repre- 
sentation of a single straight waveguide for a fixed propa- 
gation length of 40 /im. There is a perfect agreement be- 
tween the numerically calculated Fourier transform and 
its analytical solution given by Eq. [9l For integration 
boundaries greater than the propagation length the two 
calculations match with the Lorentz function generally 
used. 

Since the momentum space representation of a straight 
waveguide can be expressed for any integration bound- 
aries, we now look for the expression of a curved waveg- 
uide alone. On Fig.[8](d), the — 7r/4 profile shows the ex- 
perimental momentum space of a 90° curved waveguide 
linked by an input and an output waveguide. We ap- 
proximate in the following a solution considering a fully 
symmetric configuration consisting of a lossless circular 
waveguide. Losses are omitted to maintain the radial 
symmetry valid between [0;27r]. The Fourier transform 
in a 360° bend can be expressed in polar coordinates as 
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FIG. 9. Effect of the DLSPPW length on Fourier plane cal- 
culations. Dashed lines: calculated Fourier plane (FT) cross 
section for L = 0.5 x L spp and L = 10 x L spp . Crosses: an- 
alytical solutions (AS) calculated from Eq. [9] with the same 
integration boundaries. Red solid line: calculated Lorentzian 
profile of the wave- vector distribution of a plasmon mode with 
a propagation length of 40 /im. When L < L spp Gibbs oscil- 
lations and widening of the wave- vector content occur. 



H r (r,9) is given by Eq. |4j Since this function is pe- 
riodic with respect to between [0;27r], we can expand 
H r (r,6) in a Fourier series: 



+oo 



H r (r,6)= H n(r)e ine , 



n= — oo 



where the n th harmonic is 



1 



H n {r) = j- I H r (r,0)e- mt, d0 



(11) 



(12) 



-Mc) 




sinc(fcb.R c 7r — rar). 



The Fourier transform of the field (Eq. [T0|) can now be 
written as follow 



H n {k r ,4>)= V H)V* n 27r / rH n (r)J n {2nk r r)dr, 

Jo 

(13) 



n= — oo 



with 



and 



H n (k r ,(j)) — * h n f n (k r 



h n = smc(kbR c 7T — rar), 



(14) 



(15) 



(10) 



f n (k r ,<j>) = (-i) n e^ n 2n re 
Jo 



-(r-Rc) 



■ J n (2iTk r r)dr. 

(16) 
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FIG. 10. (a-d) Selected region of interest of the experi- 
mental momentum-space images for curved waveguides with 
R c = 19 /im, 13 jLtm, 9 /im, and 7 /mi, respectively, (e) to 
(h) Experimental profiles (shaded areas) of the wave-vector 
distribution taken along the dashed line in (a) to (d). The 
solid blue lines are the calculated profiles using Eq. [6] already 
reported in Fig. [8] The red dashed lines show the analyti- 
cal solutions obtained with the n th order rendering the best 
agreement with the data. n=100, 68, 43, and 37, respectively. 

J n (27rk r r) is the Bessel function of the first kind of n th 
order. Further simplification can be achieved by extract- 
ing the main n th harmonic h n of the series. It is deduced 
from Eq. [15] and is such that 

kbR c 7r — mr = <=> n = ki>R c . (17) 

Therefore the Fourier series is dominated by the har- 
monic integer that is the closest to the product ki>R c . 

Figure [10] illustrates a comparison between the ex- 
perimental data extracted from momentum-space spec- 
troscopy, the Fourier calculations derived from Eq. [6] al- 



ready used in Figure [5J and the analytical approxima- 
tion discussed above. The experimental Fourier planes 
of Figs. [10] (a) to (d) only show the region of interest 
for four different radii. The corresponding wave-vector 
profiles extracted and calculated at — 7r/4 (dashed lines 
in figure [7(d)) are plotted in the graphs of Figs. [10] (e) 
to (h). For radius R c > Ri, the Fourier calculations and 
the analytical solutions are in good agreement with the 
experimental data. For large radii the signature of the 
bend in momentum space is defined only by the phase dif- 
ference kb(27rR c ) = n27r. This corresponds to resonance 
condition of a ring resonator. 

Experimentally, the periodic condition on is not re- 
spected since the structure considered is formed by an 
arc of circle. Nonetheless, using the n = 100 th order of 
the Bessel function in Eq. [13]for R c =19 /im and the 68 th 
order for R c =13 /im the experimental data and the calcu- 
lated Fourier transform can be well reproduced (Figs. [10] 
(e) and (f)). For radius R c < Ri , the calculated profiles 
deduced from Eq. [6] and from the n = 47 th and n = 37 th 
orders for R c =9 /im and R c =7 /im, respectively are devi- 
ating from the experimental cross-cuts (Figs. [10] (g) and 
(h)). This disagreement is expected since none of the two 
models (Fourier transform and analytical approximation) 
are including bending loss. 

CONCLUSION 

By using dual-plane leakage radiation microscopy 
we have fully quantified the key parameters charac- 
terizing two important dielectric-loaded surface plas- 
mon polariton routing devices: linear couplers and 90° 
curved waveguides. We unambiguously demonstrated 
the added-value of performing momentum- space spec- 
troscopy. The degeneracy lift for strongly coupled waveg- 
uides and the symmetry of the split modes can be directly 
visualized and quantified. The wave- vector distribution 
associated to the curved section of the waveguide was 
also readily observed. We developed a numerical and an 
analytical analyzis to understand the experimental mo- 
mentum distribution. We found that for large radii (van- 
ishing bending loss), we can link the plasmon signature 
in Fourier space with the geometrical and modal proper- 
ties of the bend structure. The radial dependence of the 
wave- vector distribution is governed by the phase differ- 
ence kb(27rR c ). For smaller radii of curvature the bend 
loss need to be accounted for by developing an approach 
including realistic field shape. 
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